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CO ■ 1 Abstract 



The effect of the zero centered Gaussian random magnetic field distribution 
on the phase diagrams and ground state magnetizations of the transverse Ising 
thin film has been investigated. As a formulation, the differential operator tech- 
nique and decoupling approximation within the effective field theory has been 
used. The variation of the phase diagrams with the Gaussian distribution width 
(a) has been obtained and particular attention has been paid on the evolution 
of the special point coordinate with distribution parameter. In addition, the 
ground state longitudinal and transverse magnetization behaviors have been 
investigated in detail. 

2 Introduction 

Recently there has been growing interest both theoretically and experimentally 
on the finite magnetic materials especially on semi-infinite systems and thin 
films, since the magnetic properties of free surfaces are drastically different from 
the bulk material, because the free surface breaks the translational symmetry, 
i.e. surface atoms are embedded in an environment of lower symmetry than 
that of the inner atoms [U [2] . If the surface exchange coupling is greater than 
a critical value, the surface region can exhibit an ordered phase even if the bulk 
is paramagnetic and it has a transition temperature higher than the bulk one. 
This fact has been observed experimentally O 01 [5]. 

Also, the development of the molecular beam epitaxy technique and its appli- 
cation to the growth of thin metallic films has stimulated renewed interest in thin 
film magnetism. It was experimentally found that, the Curie temperature and 
the average magnetic moment per atom increases with the increasing thickness 
of the film [HE]. Thin films can be modeled by Ising model due to the fact that, 
many ultrathin films do indeed exhibit a strong uniaxial anisotropyj8]. These 
systems have been widely studied in literature by means of several theoretical 
methods such as Monte Carlo (MC) simulations [9j [TOJ [TTJ QH 031 03] , mean 
field approximation (MFA) [HI [16] and effective field theory (EFT) (171 [18]. 

But the study of ferroelectric films treated by transverse Ising model (TIM) 
are more common in the literature. Oxide thin films with pcrovskitc-type struc- 
ture (such as Barium titanate BaTiOj, and Strontium titanate SrTiO^) have 
been fabricated [TH1 US HI] and novel physical properties found in these ferro- 
electric thin films which have many application areas in technology |22j , such as 
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optoelectronics [53] , microelectronics [53] and random access memory [55] . The- 
oretically, the TIM has often been applied for the investigation of these systems. 
TIM was originally introduced by Gennes et al. [55J to describe the bulk phase 
transition of hydrogen-bonded ferroclectrics such as KH2PO4. Since then, this 
model has been applied to several physical systems like DyV02 ■ After then, this 
model has been extended to study of several other systems, such as JahnTcllcr 
and ferromagnetic systems [27] treated by means of a wide variety of techniques 
including MFA [53 [5J [301 133 131 133 133] and EFT [35] [3J [33 [33 132] • Effect of 
the dilution on the phase diagrams and magnetic properties of the ferroelectric 
thin films (by means of the TIM) has been studied within the EFT formulation 
[33 [3D 133 S3 E3] ■ Several authors have also taken into account the long range 
interactions [33E3E2]- In order to simulate more realistic models, different ex- 
change interactions have been defined and different transverse fields have been 
considered on the surface and bulk, which anticipates the mimics of the surface 
effects within the framework of EFT [33[33[53[5I][S3[53[M1[S5]- Another 
improvement of the Ising thin film models is the consideration of amorphisation 
of the surface due to environmental effects, and these realizations have also been 
solved with using EFT [53 [57] [5_3 E3 [63 EI] . Moreover, in order to mimic the 
effect of some imperfections on the film characteristics, straight line of magnetic 
adatoms [52], defect layers [53] and seeding layer effects [53] have been stud- 
ied. Apart from these works, decorated transverse Ising thin films with EFT 
[53 [53 [57] [55] and MC [5j|] have been studied. There are also higher spin Ising 
thin films e.g. spin-1 Ising thin films have been studied [70l [7T] . 

In the semi-infinite systems, depending on the ratio between surface ex- 
change interaction and bulk exchange interaction, the system may order on the 
surface before it orders in the bulk which is called extraordinary transition. In 
the contrary of this, ordinary transition means that the surface critical temper- 
ature is the same as the bulk transition temperature. In the phase diagram, 
the intersection point between these two transitions is called special point. As 
the film gets thicker, it approaches the semi-infinite system; this unusual effect 
shows itself as an intersection point in the phase diagrams of the films that 
have different thickness plotted in critical temperature versus surface exchange 
interaction plane. 

On the other hand, the Ising model in a quenched random field (RFIM) has 
been studied over three decades. The model which is actually based on the 
local fields acting on the lattice sites which are taken to be random according 
to a given probability distribution was introduced for the first time by Larkin 
[75] for superconductors and later generalized by Imry and Ma [73] ■ Beside the 
similarities between diluted antifcrromagncts in a homogenous magnetic field 
and ferromagnetic systems in the presence of random fields [73] [75] , a rich class 
of experimentally accessible disordered systems can be described by RFIM such 
as structural phase transitions in random alloys, commensurate charge- density- 
wave systems with impurity pinning, binary fluid mixtures in random porous 
media, and the melting of intercalates in layered compounds such as Tz&*[76]. 
RFIM generally models the phase transitions and interfaces in random media 
|77[ 178] , e.g prewetting transition on a disordered substrate can be mapped to 
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2D RFIM problem[7g. 

The aim of this work is the determine the Gaussian random longitudinal 
magnetic field distribution on the phase diagrams and magnetization behavior of 
the transverse Ising thin film. For this aim, the paper is organized as follows: In 
Sec. |3]we briefly present the model and formulation. The results and discussions 
are presented in Sec. El and finally Sec. [5] contains our conclusions. 

3 Model and Formulation 

Thin film can be modeled by a layered structure which consist of interacting 
L parallel layers. Each layer is defined as a regular lattice with coordination 
number z. The Hamiltonian of the thin film is given by 

w = - E J ^ s i "E s > - E^ s * C 1 ) 

<i,j> i i 

where sf is the z component of the spin-i operator (which can take the values 
sf ± 1) and sf is the x component of the spin-| operator at a lattice site i. The 
transverse field Q is assumed to be distributed regularly on the lattice sites while 
the local external longitudinal magnetic field Hi on a lattice site i, is distributed 
via a given probability distribution. The exchange interaction Jij between the 
spins on the sites i and j takes the values according to the positions of the nearest 
neighbor spins. The two surfaces of the film have the intralayer coupling J\. 
The intcrlayer coupling between the surface and its adjacent layer (i.e. layers 
1, 2 and L — 1,L) is denoted by J 2 . For the rest of the layers, the interlayer and 
the intralayer couplings are assumed as J 3 . The first summation in Eq. ([1]) is 
over the nearest-neighbor pairs of spins, and the remaining summations are over 
all the lattice sites. Magnetic fields are distributed on the lattice sites according 
to a given Gaussian probability distribution 

where a is the width of the zero magnetic field centered Gaussian distribution. 

Number of L different representative longitudinal magnetizations in z di- 
rection (to, = (sf ) , i = 1, 2, . . . , L) for the system can be given by usual EFT 
equations which are obtained by differential operator technique and decoupling 
approximation (DA) [501 1ST] . 

mi = [Ax + miBiY [A 2 + m 2 B 2 ] 

m L = [Ax + m L B{\ z [A 2 + TOi_iS 2 ] 

m 2 = [As + m 2 B 3 ] z [A 2 + m Y B 2 \ [A 3 + m 3 B 3 ] 

m L -i = [A 3 + m L -iB 3 ] z [A 2 + m L B 2 ] [A 3 + m L ^ 2 B 3 ] 

m k = [A 3 + m k B 3 ] z [A 3 +m k ^ 1 B 3 ][A 3 +m k+1 B 3 },k = 3,4:,...,L-2. 

(3) 
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Here m,, (i = 1,2, ... ,z) denotes the longitudinal magnetization of the i th layer. 
The coefficients in the expanded form of Eq. ([3]) are given by 



jdH t P (Hi) cosh* (J P V) cosh' (J 3 V) sinh m (J p V) sinh" (J,V) / (fli.x) | x=0 



(4) 

where V is the usual differential operator in the differential operator technique 
with indices p,q = 1,2,3. The function / (Hi,x) is defined by 



f(H h x) 



x + H.; 



tanh(/3y), y= (x + H t j 2 + ti 2 



1/2 



(5) 



as usual for the spin-i system with transverse field. In Eq. ([5]), (3 = l/(/csT) 
where ks is Boltzmann constant and T is the temperature. The effect of the 
exponential differential operator to an arbitrary function F(x) is given by 



exp (aV)F (x) = F (x + a) 



(6) 



with any constant a. DA will give the results of the Zernike approximation [82) 
for this system. Transverse magnetization for the i th layer (mf = (sf) , i 
1.2 I . . can be obtained by using the function 



f x (H t ,x) = - tanh (py) , y=Ux + H t ) 2 + Q 2 

y L 



1/2 



(7) 



in Eq. @ (and then in Eq. ©) instead of Eq. ([5]). 

With the help of the Binomial expansion, Eq. ([3]) can be written in the form 

z 1 

mi = ££#1 m i m 2 



i=0j=0 

z 1 



m L 



m 2 



m L -i 



2 1 1 



i=0 j=0fc=0 
z 1 1 



= EES** (»'»**) 



m L-l m L m L-2 



(8) 



i=0j=0fc=0 

Z 1 1 

EEE^ 3 (*> j ' fc ) m fe m i-i m fe+i 

i=0 j=0fc=0 



where 



#3 (*, j, k) 



A l f 1 A z + 1 - l - k B 3 2 Bl +k 



(9) 
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These coefficients can be calculated from the definitions given in Eq. (01 with 
using Eqs. ([5]) and J6|). 

For a given Hamiltonian and field distribution parameters, by determining 
the coefficients from Eq. Q we can obtain a system of coupled non linear 
equations from Eq. ([8]). and by solving this system we can get the longitudinal 
magnetizations of each layer (m,,i = 1,2, ...,L). The longitudinal magneti- 
zation (m) and transverse magnetization (m x )oi the system can be calculated 
via 

1 L 1 L 

m = -j^Tm l , m x = j^J 71 *' (10) 

i=i i=i 

Transverse magnetization of each layer (mf ) can be obtained by constructing 
non linear equation system given in Eq. ([8]) by calculating coefficients from Eqs. 
jH and ([9]) using Eq. ([7]), instead of Eq. ©. 

Since all longitudinal magnetizations are close to zero in the vicinity of the 
critical point, we can obtain another coupled equation system to determine the 
transition temperature by linearizing the equation system given in Eq. ([8]), i.e. 

mi = Ki (1,0) mi + Ki (0,1) m 2 

m L = K 1 (l,0)m L + K 1 (0,l)m L -i 

m 2 = K 2 (1, 0, 0)m 2 + K 2 (0,1,0) mi + K 2 (0,0,l)m 3 (11) 

m L _i = K 2 (1, 0, 0) m L - 1 + K 2 (0,1,0) m L + K 2 (0,0,1) m L - 2 

m k = K 3 (l,0,0)m k + K3(0,l,0)m k - 1 +K 3 (0,0,l)m k+1 . 

Critical temperature, as well as critical transverse field can be determined 
from dct(A) = where A is the matrix of coefficients of the linear equation 
system given in Eq. ([TT]) . 



4 Results and Discussion 

In this section we discuss the effect of the Gaussian random magnetic field 
distribution on the phase diagrams of the system. We restrict ourselves in the 
case that all interactions are ferromagnetic i.e. Ji > 0, (i = 1,2,3) and z = 4 
i.e. all layers in thin film have square lattices. We use the scaled interactions as 

Js = J, r n = y , n = 1, 2 n = 1 + A s (12) 

Note that, the parameter A s mostly used in literature in order to investigate 
the system according to a surface exchange interaction in relation with other 
exchange interactions in the system. Since we want to concentrate on the effect 
of a and A s on the phase diagrams and magnetizations, we choose r 2 = 1.0 in 
this work. 

4.1 Phase Diagrams 

Gaussian distribution defined in Eq.(f2]) is governed by only one parameter a , 
which is the width of the distribution. This distribution distributes negative 
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and positive valued magnetic fields -which are chosen from the zero centered 
Gaussian distribution- to lattice sites, so that the overall magnetic field is equal 
to zero. Although the total magnetic field is zero, randomly distributed nega- 
tive and positive fields drag the system to the disordered phase. On the other 
hand, the interactions J„, (n = 1, 2, 3) enforce the system to stay in the ordered 
phase. Another factor is the transverse field which forces the spins to align 
perpendicular to z direction. The last factor, the temperature which causes 
thermal agitations induces a disordered phase when the energy supplied by the 
temperature to the system is high enough. Thus a competition takes place 
between these factors. 

The effect of the Gaussian magnetic field distribution on the phase diagrams 
of the bulk systems is now well known. As a rises, the critical temperature 
decreases continuously which originates from the rising randomness of the ori- 
entations of the s| . Similarly, rising 51 forces the spins to align perpendicular to 
the z direction, and consequently it decreases the critical temperature. Transi- 
tion temperature reduces to zero at the critical value of transverse field f2 = fi c . 
After this value of f2 c , the system cannot exhibit an ordered phase, due to the 
lack of the chance of the spins to front in z direction. 

On the other hand, the phase diagrams for the Ising model with thin film 
geometry in the (fcsT c /J, A s ) and (f2 c /J, A s ) plane are also well known (see 
e.g. Ref. [HI]). A s ^ s rises, the critical temperature and the critical transverse 
field of the system increases for arbitrary L. The curves for the different film 
thickness (L) intersect at a special point which can be denoted by (A*, ksT* / J) 
and by (A*,fi*/J) for the curves in the (fcsT c / J, A s ) and (Cl c /J,A S ) planes, 
respectively. For the values of A s > A*, thin films have higher critical temper- 
ature (and critical transverse field) than the bulk system while the reverse is 
valid for A s < A*. For the value of A s = A*, thin film with thickness L has 
the same critical values of the bulk, independent of the L. The value of A s also 
changes the relation between the film thickness and the critical temperature 
(and critical transverse field) of the film. The thicker films have higher critical 
values for A s < A* while thicker films have lower critical values for A s > A* 
than thinner ones. 

After this short summary, at first let us investigate the effect of the Gaussian 
random field distribution width on the phase diagrams of the transverse Ising 
thin film in (fcsTc/J, A s ) and (fi c /J, A s ) planes. In Figs. ([1]) and ([2]), phase 
diagrams can be seen for different L values with chosen a. We can see from 
Figs. ([1]) and ([2]) that, rising a decreases the critical temperatures and critical 
transverse fields as in the bulk counterpart of discussed system. This decrement 
in the values of the UbTc/ J and £l c / J shows itself also in the special point 
coordinates ksT* / J and 0*/J while there is no significant change in the value 
of the A* when a rises. We note that for a = 0.0, A* value that makes critical 
temperature at f2 = 0.0 independent of the film thickness (L) is A* = 0.3068 
and the value of A* that makes critical transverse field independent of the film 
thickness (L) is A* = 0.3326 as in Ref. [49] . Also, rising a does not affect the 
relation between the film thickness and critical temperature (or transverse field) 
, i.e. for A s < A* thicker films have higher critical values and reverse is true for 
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A s > A* with a ^ 0, if the special point is present. When a is large enough, 
special point disappears and thicker films have lower critical values for a certain 
A s as seen in Figs, ©(f) and ©(f). 

The variation of the ksT* / J and SI*/ J with a can be seen in Fig. ©. 
This variation is nothing but the variation of the ksT c / J (at SI = 0) and Sl c / J 
with a for the TIM in the bulk system, since the special point is defined as the 
intersection point of the curves in the (/cbT c / J, A s ) (and (Sl c /J, A s )) plane for 
different film thicknesses and with the bulk curve in the same plane, which is 
just a parallel line to the A s axis. We note that the curves in Fig. © ends at 
the point a = 4.22, i.e. for the values that provide a > 4.22, the system can not 
exhibit special point. 

The phase diagrams in (ksTc/ J, SI/ J) plane can be seen for the values of 
A s = 0.0 and A s = 0.6 (which provide the conditions A s < A* and A s > A* 
respectively) with film thicknesses L = 4, 6, 10 and for different a values in Fig. 
©. We can see from Fig. © that the qualitative relation between critical 
temperatures (and critical transverse fields) and film thickness is also valid for 
these diagrams; i.e. for A s = 0.0 < A*, all critical temperatures increase 
whereas for A s = 0.6 > A*, all critical temperatures decrease when the film 
thickness increases for given A s and any a . As in the bulk systems, rising a 
shrinks the ferromagnetic region in the (ksT c / J, SI/ J) plane. This fact is due to 
two reasons: rising SI forces the spins to align perpendicular to the z direction 
and rising a forces the spins to align parallel to z direction, but randomly i.e. +z 
or —z. Two of these results to the decline of the order parameter. Now, in order 
to understand how a effects the phase transition beside the transverse field, let 
us look at the effect of these parameters on the longitudinal and transverse 
magnetization. 

4.2 Magnetization Behaviors 

Let us investigate the effect of the random Gaussian magnetic field distribution 
on the magnetization behaviors with temperature for both the longitudinal and 
transverse counterparts. In order to make inferences about how a effects the 
variation of longitudinal magnetization with temperature, we present the curves 
for the parameter values L = 4, 10, SI = 0.0, 2.0 and A s = 0.0, 0.6 through large 
set of parameters as representatives. These curves can be seen in Fig. © 
for the film thickness L = 4, and in Fig. © for the film thickness L = 10. 
Each plot contains surface longitudinal magnetization (mi) (with solid lines) 
and inner layer longitudinal magnetization (rriL/2) (with dotted lines), in order 
to investigate how this randomness spreads into the film. 

Firstly, there is no need to elaborate on the general effect that, rising a 
decreases the longitudinal magnetization at any temperature that lower than 
the critical temperature. At first glance, from Figs. ©(a),(c) and © (a),(c), we 
can see that, in the absence of the transverse field, rising a separates the ground 
state longitudinal magnetizations of the surface and inner layer, for a given L 
and A g . The value of the difference between the ground state longitudinal 
magnetization of the surface and the inner layer depends on a, L and A s in 
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the absence of the transverse field. As we can see from Figs, ©(b), (d) and 
([6]) (b),(d), the same effect can be achieved by rising CI even in the absence 
of random field distribution (see the curves labeled by A in Figs. (0 (b), (d) 
and (0 (b), (d)). Rising CI results in that more spins align perpendicular to z 
direction. In the case of A s = 0.0, i.e. all exchange interactions are equal to 
each other, since the spins which are located in the inner layers have one more 
neighbor spin than the spins located at the surface layer; it is expected that the 
inner spins have more chance to align in z direction than that of the surface 
spins. Accordingly, we get higher longitudinal magnetization in the inner layers 
than that of the surface layer (the curve labeled by A in Figs. ([5]) (b) and (0 
(b)). Rising of A s compensates the lack of a neighbor site of the surface spins, 
then after some value of A s , the surface layer magnetization becomes greater 
than that of the inner ones (the curve labeled A in Figs. (0 (d) and ([6]) (d)). 
The same scenario is also valid for the random field case. Rising a causes the 
random alignment of all spins in ±z direction in the system. Because of the 
existence of an additional spin interaction, the spins which are located in the 
inner layers have more chance to align parallel to each other than the spins 
in the surface layer. When A s increases, after some value of it, this situation 
turns reverse in the presence of Gaussian random field distribution (for instance, 
compare the curves labeled C of Figs. (0 (a) and ([5]) (c)). 

Increasing a causes an increment in the difference between the longitudinal 
magnetization of the surface layer and the inner layer at any temperature (that 
little than the critical temperature) and transverse field value (see for example 
the curves labeled C and D in Fig. (0(c)). As L increases, this difference also 
increases for large enough a values (see for example the curves labeled D in Fig. 
(0(c) and ©(c)), and rising A s value at first reduces this difference and after a 
specific value at which the longitudinal magnetization of the surface and inner 
layer become equal to each other, rising A s increases the difference between the 
two magnetization profiles, and the longitudinal magnetization of the surface 
becomes greater than that of the inner one. 

In order to focus on the difference between the ground state longitudinal 
magnetization of the surface (mi) and that of the inner layer (m^ 2 )i we present 
the variation of the Am = mi — m L / 2 curves with the transverse field at the 
temperature ksT c /J = 0.001 for some selected values of A s , L and a in Fig. ([7]). 
At this temperature, the energy supplied from the temperature to the system is 
very small in comparison with the spin-spin interaction, then the situation may 
be treated as the investigation of ground state magnetization. In the absence 
of the randomly distributed field (i.e. a = 0.0), Am ^ value originates 
completely from the transverse field. The general trend of variation of Am with 
CI can be seen in Figs. (0 (a), (e) in the absence of the magnetic field. Rising 
Cl can not create significant change until Cl/J = 1.0, then |Am| increases with 
increasing Cl/J until the specific value of Cl/J for a given A s and L. After 
then, |Am| approaches to zero and finally at the value of Cl c /J for a given A s 
and L, |Am| falls to zero because all layers exhibit zero magnetization. We can 
see from Fig. ([7]) that this trend is also valid in the presence of the randomly 
distributed longitudinal magnetic field. For instance, we can see from Fig. (0 



(b) that the shapes of the curves are similar to that given in Fig. (J7]) (a) but 
with lower Cl c / J and higher value of | Am| at any Cl/ J (even in the absence of the 
transverse field) for fixed A s . This fact comes from the randomly distributed 
field which has mechanism explained above. But undoubtedly, the net effect 
of the transverse field and randomly distributed magnetic field on the Am is 
not simply the sum of them, since the random magnetic field distribution also 
affects the transverse magnetization which originates from the transverse field. 
We can make one more conclusion that, the higher a curves have slightly higher 
maximum value of the |Am|, which can be seen in Figs. (|7|) (c) and (d) by 
comparing the curves labeled by A,B,C,D. 

Finally, we want to make some conclusions on the effect of the a on the 
ground state transverse magnetization of the surface layer and an inner layer. 
For this purpose, let us investigate the situation for the values of Cl that provides 
the condition Cl/J > 1, since only these values makes Am different from zero 
in the absence of the magnetic field, as seen in Fig. (J7J (a). We again focus 
on the Am 1 = mf — TO £/2 calculated at temperature ksT c /J = 0.001. We 
can see the variation of the mf and m x h , 2 with a for L = 4 and Cl = 2.0 and 
A s = 0.0,0.4 in Fig. ((U) (a). We can see that, rising a decreases both mf and 
m L/2 va hics, which can be explained by the fact that when a rises, more and 
more spins are forced to align in ±z direction and this causes a decline in the 
transverse magnetization. As discussed above, Am x is positive at A s = while 
it turns to negative when A s = 0.4 at any a, which can be seen in Fig. ([5]) 
(a). If we look at the variation of the mf — a and ™^/2 — a curves calculated 
at a temperature fcsT c /J = 0.001 for Cl = 1.3 and L = 4, we can see from 
Fig. ([5]) (b) that some mf curves (for the values A s = 0.14,0.16) are located 
above the m x L , 2 curve , and some of them (for the values A s = 0.24, 0.26) are 
below the m|y 2 curve f° r au values of a. But the curves of mf for the values 
A s = 0.18,0.20,0.22 intersects the curve Tnf L , 2 \ while a rises. This means that 
for some values of A s , Am x is always positive or negative (i.e. for all a), but 
for some values of A s , Am x changes its sign as a varies. These A s values are 
of course depend on the other parameters of the system such as L, Cl. The 
borders of the positive valued Am 1 and negative valued Am 1 (calculated at the 
temperature fcsT c / J = 0.001) can be seen in Fig. © for L = 4 and for selected 
values of A s = 0.18,0.20,0.22,0.24 in the (a, CI /J) plane. Ground state Am 1 
calculated at A s is positive on the region that covers the right side of the related 
A s curve in the (a, CI/ J) plane and negative on the region that covers left side 
of the related A s curve. We can see from the Fig. ^ that, rising A s shifts the 
curves to the right and this means that expanding the region that Am 1 < 0. 
The absence of A s < 0.18 curves in Fig. ^ tells us that, for A s < 0.18 always 
Am 1 > within the region in Fig. © and similarly, the absence of A s > 0.24 
curves in Fig. (J9)) means that, for A s > 0.24 always Am 1 < within the 
region in Fig. (j9j). For a fixed value of er, and for example with A s = 0.22, 
when Cl increases, the system passes from the region Am 1 < to the region 
Am 1 > at a specific value of Cl/ J. But, for instance for a fixed value of 
Cl/J = 1.6, rising a crosses from the region Am x < to the region Am 1 > 
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for a certain er, then with rising a crosses Am 1 < again, according to the 
curve related to the A s = 0.22 (which is labeled as C in Fig. ([9])). This shows 
us that, indeed a competition between the Gaussian distributed longitudinal 
magnetic field and the transverse field, results some complicated situations on 
the transverse magnetization. 

5 Conclusion 

In this work, the effect of the zero centered Gaussian random magnetic field 
distribution on the phase diagrams and ground state magnetizations of the 
thin film described by the TIM has been investigated. As a formulation, the 
differential operator technique and DA within the EFT has been used. 

First, the evolution of the phase diagrams of the thin film in the (ksT c /J, A s ), 
(fl c /J, A s ) and (fcsT c /J, 0,/J) planes with random Gaussian longitudinal mag- 
netic field distribution width (er) investigated in detail. Also the effect of the a 
on the special point which appears in the phase diagrams in the (ksT c / J, A s ) 
and (f2 c /J, A s ) planes has been investigated. For the coordinates of the special 
point (A*, k B T* / J) and (A*, Cl*/J) in the (k B T c /J, A s ), (Q c /J, A s ) planes, re- 
spectively; we found that rising a makes no significant change in the A* in both 
planes, but gives rise to a decline in the fcsT* / J and fl*/J- After the value of 
a = 4.22, the special point disappears in both planes simultaneously Although 
the special point can not appear for the values a > 4.22, thin film can achieve 
ordered phase after some value of A s as seen in Figs. [T] (f ) and[2] (f ) . This region 
shows itself in the corresponding semi-infinite system with ordered surface and 
disordered bulk. At the same time, the effect of the rising a on the relation 
between the critical temperature/transverse field and film thickness has been 
discussed in detail. 

Although there are different physical mechanisms underlying the observa- 
tions, both rising 57 and a give rise to a decline in the critical temperature. 
In order to focus on the effect of a on the phase transition in the presence of 
transverse field, the variation of the ground state longitudinal/transverse mag- 
netizations with a has been obtained. After making some conclusions on the 
effect of the rising a on the curves that represents the evolution of the magneti- 
zation with temperature, we have focused on the effect of the er on the ground 
state magnetizations of the surface and inner layer. It is not surprising that, for 
small A s , the surface ground state transverse magnetization is greater than the 
inner one in the absence of the magnetic field. As explained in the Sec. 0] rising 
A s reverses this situation, i.e. in terms of the defined parameter in the Sec. HI 
Am x = mf — m£/ 2 < which is calculated at the temperature UbTc/ J = 0.001. 
It has been shown that, within the (er, J7/J) plane, for some values of the A s , 
we get Am x > while for the higher values we always have Am x < regardless 
of er and fi value. However, in the region that covers the transition between the 
two cases, there is a complicated situation which comes from the competition 
between and er. 

We hope that the results obtained in this work may be beneficial form both 
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theoretical and experimental point of view. 
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Figure 1: The phase diagrams of the thin film described by the TIM with single 
Gaussian random field distribution in the (fcsT c / J, A s ) plane for different a 
values with selected thickness values of L = 4, 6, 10 . The selected value of r2 is 
r<2. = 1 
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Figure 2: The phase diagrams of the thin film described by the TIM with single 
Gaussian random field distribution in the (O c /J, A s ) plane for different a values 
with selected thickness values of L = 4, 6, 10 . The selected value of T2 is r2 = 1 
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Figure 3: Variation of ksT* / J and f2*/J with a for the thin film described by 
the TIM with single Gaussian random field distribution. 
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Figure 4: The phase diagrams of the thin film described by the TIM with single 
Gaussian random field distribution in the (fcsT c / J, Q,/ J) plane for different a 
values with selected thickness values of L = 4, 6, 10 . The chosen value of ri is 
r% = 1 
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Figure 5: Variation of the longitudinal layer magnetizations with temperature 
for some selected values of A s , O and a for the thin film described by the TIM 
with single Gaussian random field distribution with thickness L — 4. Surface 
longitudinal magnetization m\ is represented by solid line while the longitudinal 
magnetization of the L/2 indexed layer (m L / 2 ) is represented by dotted line. 
The chosen value of r% is r2 = 1. 
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Figure 6: Variation of the longitudinal layer magnetizations with temperature 
for some selected values of A s , Q and a for the thin film described by the TIM 
with single Gaussian random field distribution with thickness L — 10. Surface 
longitudinal magnetization mi is represented by solid line while longitudinal 
magnetization of the L/2 indexed layer (m L / 2 ) is represented by dotted line. 
Chosen value of r 2 is r 2 = 1 . 
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Figure 7: Variation of the difference between the longitudinal magnetization 
of the surface layer and the inner layer (Am = mi — m L / 2 ) calculated at the 
temperature &bT c /J = 0.001 with the transverse field Vt/J for some selected 
values of A s , i and a for the thin film described by the TIM with single Gaussian 
random field distribution. The chosen value of r 2 is r 2 = 1. 
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Figure 8: (a): Variation of the surface layer transverse magnetization (mf, 
solid line) and inner layer transverse magnetization (w^y 2 , dotted line) with 
a evaluated at the temperature fcgT c / J = 0.001 for the values of A s = 0.0 
and A s = 0.4. Fixed values are r-i = l,L = 4 and il = 2.0. (b): Variation 
of the surface layer transverse magnetization {mf, dotted line) and inner layer 
transverse magnetization (?n|, 2 , solid line) with a evaluated at the temperature 
k B T c /J = 0.001 for the values of from A s = 0.14 to A s = 0.26. The increment of 
A s from one curve to the neighbor curve is 0.02. Fixed values are r 2 = 1,1 = 4 
and fl = 1.3. 
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Figure 9: The borders of the positive valued Am 1 and negative valued Am x 
calculated at the temperature fcsT c /J = 0.001 for the film thickness L = 4 and 
for selected values of A s = 0.18,0.20,0.22,0.24 in the (a.Q/J) plane. Am x is 
positive (i.e. mf > rn x L / 2 ) 011 tne re gi° n that covers the right side of the related 
A s curve and negative (i.e. mf < m x L ^ 2 ) on the region that covers the left side 
of the related A s curve. 
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